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We investigate a Hamiltonian lattice version of the two-dimensional Wess-Zumino model, with special emphasis 
to the pattern of supersymmetry breaking. Results are obtained by Quantum Monte Carlo simulations and Density 
Matrix Renormalization Group techniques. 



1. INTRODUCTION AND PRESENTA- 
TION OF THE MODEL 

Most studies of lattice field theory are per- 
formed in the well-known Lagrangian formalism, 
discretizing both space and time. We feel that it 
is important not to neglect the Hamiltonian for- 
malism which affords, e.g., a more immediate 
approach to the mass spectrum and to the im- 
plementation of fermions; numerical methods are 
quite different in the Lagrangian and in Hamilto- 
nian formalism, and comparison between the two 
can provide important test of the methods as well 
as of universality. 

We chose to study an Hamiltonian lattice ver- 
sion of the two-dimensional Wess-Zumino model, 
described in Ref. [2]; the model enjoys an exact 
1-dimensional supersymmetry algebra H = Q^, 
which suffices to preserve a few key properties 
of the continuum model (which has a full 2- 
dimensional supersymmetry): the energy spec- 
trum is non- negative; states of nonzero energy ap- 
pear in boson- fermion doublets; spontaneous su- 
persymmetry breaking is equivalent to a strictly 
positive ground-state energy Eq. 

The model is parametrized by the prepotential , 
an arbitrary polynomial in the bosonic field V{(l)); 
the model is superrenormalizable, and the only 
renormalization needed is the normal ordering of 
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It is interesting to notice that the strong- 
coupling expansion predicts supersymmetry 
breaking if and only if the degree of the pre- 
potential V is even, while perturbation theory 
predicts supersymmetry breaking if and only if 
V has no zeroes; the two contitions are quite 
different in the case, e.g., of quadratic V, and it 
is interesting to study the crossover from strong 
to weak coupling. 



2. NUMERICAL METHODS 

A standard numerical simulation technique 
that can be applied to the model we are interested 
in is the Green Function Monte Carlo (GFMC) 
algorithm ^^i,- The interested reader may find a 
description of GFMC applied to the present prob- 
lem in Ref. '7; we only remark here that the 
GFMC algorithm generates a stochastic repre- 
sentation of the ground-state wavefunction, from 
which expectation values of observables can be 
computed; the variance of observables can be kept 
under control by the introduction of a population 
of K walkers (field configurations), and an ex- 
trapolation to K —>■ oo oi the result is required. 

Another numerical approach which has been 
applied with success to 1-dimensional Hamilto- 
nian models is the Density-Matrix Renormaliza- 
tion Group (DMRG) |4I5| : the parameters of the 
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algorithm are the dimension S to which the lo- 
cal (1-site) Hilbert space is truncated, and the 
number AI of eigenstates of the density matrix of 
the full system that are retained; extrapolation 
to M ^ oo and 5* ^ cxd is needed. DMRG is a 
deterministic algorithm and does not suffer from 
the sign problem; it is possible to obtain directly 
infinite- volume results. 

3. NUMERICAL RESULTS 

The most interesting property of the model 
is the pattern of supersymmetry breaking, es- 
pecially in the case of quadratic prepotential; 
we study the case of even prepotential V{(f>) = 
A20^ -I- Ao, for which the strong-couphng expan- 
sion predicts supersymmetry breaking for all val- 
ues of Ao, while perturbation theory predicts un- 
broken supersymmetry for Aq < 0. 

The model enjoys an approximate (pm — > ~4>rrn 
n-m — + 1 — rim symmetry ( "parity" ) ; supersymme- 
try breaking is incompatible with parity break- 
ing [H], and it is reasonable to expect only two 
phases: unbroken supersymmetry and broken 
parity for Aq < Ac, broken supersymmetry and 
unbroken parity for Aq > Ac. 

An analysis of the approach to the continuum 
limit, along a trajectory of "constant physics" 
in the phase of broken supersymmetry was pre- 
sented in Ref. In the present paper, we study 
the model varying Aq at fixed A2 — 0.5. 

We monitor supersymmetry breaking by look- 
ing at the ground-state energy density Eq/L. 
Infinite-volume DMRG results are shown in 
Fig. n we remind that they must be extrapolated 
to 5 ^ 00 and M ^ 00. GFMC results are ex- 
trapolated to L ^ cxD and if — > c» by a best fit 
to the form 

Eq a b c d 

which gives a good (at least for Aq > —0.46); 
the fit gives v ~ 0.5, signaling a slow convergence 
to iiT 00 (to stabilize the fit, we ran with K 
up to 5000 for our smallest lattice L — 22). £ is 
plotted in Fig. 12 together with DMRG results for 
the highest values of S and M; the agreement is 
remarkable. 



■ ■ S = 8. M = 10. 15, 20 

♦ ♦ S = iO, A/ = 15, 20. 25 
A A S = 12, M = 15 
< S = 16, 16 

r S= 18, M= 18 



= 0.5 




Figure 1. Ground-state energy density vs. Ao ob- 
tained by DMRG with different values of S and 
M. 
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Figure 2. Ground-state energy density vs. Ao ob- 
tained by DMRG with different values of S and 
M. 



We study parity symmetry breaking looking at 
the Binder cumulant 



B 



M : 



where the sum excludes sites closer to the border 
than (typically) 6; a good estimate of the transi- 
tion point is the intersection of the curves B vs. Ao 
obtained at different values of L. GFMC results 
are shown in Figs. |31 and ^ they are quite sensi- 
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Figure 3. The Binder cumulant B ior K = 200. 
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Figure 5. The long-distance correlation of (jj vs. 
Ao for i 58 and = 500. 
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Figure 4. The Binder cumulant B for K = 500. 



tive to K and not very precise. Therefore, we 
follow a different strategy and consider the con- 
nected correlation function Gd — {4'n4>m)c aver- 
aged over all n, m pairs with \m — n\ — d, exclud- 
ing pairs for which m or n is closer to the border 
than (typically) 6; since we are using staggered 
fermions, even and odd d may correspond to dif- 
ferent physical channels. The correlation function 
is quite different in the two phases: for Aq ^ Ac, 
there is a marked difference between even and odd 
distances, and Gd is short-ranged; for Aq Ac, 
even and odd distances are equivalent, and Gd 
appears to be long-ranged. We fit Gd to the form 



exp(a + bd) + c, separately for even and odd d; the 
best fits give a good if we remove the small- 
est distances, typically d < 3 for the odd channel 
and d < 6 for the even channel. The difference 
between the two phases is apparent, e.g., in the 
plot of c vs. Aq, presented in Fig.|Sl 
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